function res=thetapuntoe(t,X)
    global b;
    global I;
    global Kd;
    %res = zeros(1, length(X));
    res(2)=(-b/I)*X(2)+(Kd/I)*(0.01-X(2));
    res(1)= X(2);
end